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Output Feedback Control of Inhomogeneous Parabolic PDEs with Point Actuation 
and Point Measurement using SOS and Semi-Separable Kernels 

Aditya Gahlawat 1 and Matthew M. Peet 2 


Abstract —In this paper we use SOS and SDP to design output 
feedback controllers for a class of one-dimensional parabolic partial 
differential equations with point measurements and point actuation. 
Our approach is based on the use of SOS to search for positive quadratic 
Lyapunov functions, controllers and observers. These Lyapunov func¬ 
tions, controllers and observers are parameterized by linear operators 
which are defined by SOS polynomials. The main result of the paper is 
the development of an improved class of observer-based controllers 
and evidence which indicates that when the system is controllable 
and observable, these methods will find a observer-based controller 
for sufficiently high polynomial degree (similar to well-known results 
from backstepping). 

I. INTRODUCTION 

Parabolic Partial Differential Equations (PDEs) are a class of 
system used to model processes such as diffusion, transport and 
reaction. Some examples of systems which have been modelled 
using parabolic PDEs include plasma in a tokamak [22], heat prop¬ 
agation, and spatial dynamics of population in an ecosystem [14]. 
In this paper we consider a class of inhomogeneous linear scalar 
valued Parabolic PDEs. We assume that only boundary control 
and sensing is available for the PDEs. The control is exercised 
through a Dirichlet boundary condition and a Neumann boundary 
measurement of the state is available. The goal of this article is to 
use this boundary measurement to construct a boundary controller 
which ensures that the state of the system remains bounded in the 
presence of an exogenous input (has finite Z/ 2 -gain). We refer to 
this as output feedback based boundary control. 

In order to design output feedback based boundary controllers, 
we design a Luenberger observer where the error dynamics have 
finite Z /2 gain from disturbance to error. We then design a controller 
for the system which utilizes the state of the observer and show 
that for the resulting closed-loop system, there is a bound on the 
1/2 gain from disturbance to output. Our approach is based on 
parameterizing the set of quadratic Lyapunov functions by the set 
positive operators, which in turn is parameterized by polynomials 
and ultimately by Sum-of-Squares (SOS) polynomials and positive 
matrices - leading to a Linear Matrix Inequality (LMI). The 
approach we take in this paper is akin to LMI methods for control 
of linear Ordinary Differential Equations (ODEs) using Lyapunov 
inequalities and a variable substitution trick. However, because our 
inequalities are expressed as operators in Hilbert space, we refer to 
our approach as a Linear Operator Inequality (LOI). 

This article extends our work in [7] wherein we designed output 
feedback boundary controllers for a one-dimensional homogenous 
heat equation by considering a simpler class of positive operators 
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(also parameterized by SOS polynomials). This paper improves on 
the work in [7] by a) Considering the larger class of inhomo¬ 
geneous, possibly unstable parabolic PDEs b) By considering a 
larger class of Lyapunov functions defined by positive multiplier 
and integral operators with semi-separable kernels and c) providing 
evidence (but not a proof) that this new class of operators can 
be used to design output-feedback based controllers whenever the 
system is observable and controllable. Specifically, the class of 
Lyapunov functions we use has the form 

V(w) = f M(x)w(x) 2 dx+ f f w(x)K(x,^)w(^)d^dx, 
Jo Jo Jo 

where 

I<(x,0 = 

and were M, K i and K 2 are polynomials and w represents the 
spatially distributed state of the PDE. A kernel K of this form is 
referred to as semi-separable. 

One popular and relatively straightforward method for output 
feedback boundary control of PDEs is backstepping [13]. This 
method relies on constructing an invertible operator which, in 
closed loop, maps the state of the system to the state of a chosen 
stable system for which a quadratic Lyapunov function exists. Our 
approach varies in the fact that we search for both the controller 
and the quadratic Lyapunov function. Although our approach is 
different, similar to backstepping, the numerical results indicate that 
we can construct output feedback controllers for any controllable 
and observable system. Some other examples of work which use 
Lyapunov functions for analysis and control of PDEs are [4], [5], 
An example of application of LMIs for the control of PDEs is [6] 
where the authors synthesize stabilizing boundary controllers for 
uncertain semi-linear PDEs using quadratic Lyapunov functions 
parameterized by positive scalars. Early results on the use of SOS 
for analysis and control of infinite-dimensional systems can be 
found in [17], [15]. Additional recent work on the application of 
polynomials to infinite-dimensional systems can be found in the 
research done by our colleagues in [21] and [1], 

The paper is organized as follows: Section [III] outlines the 
problem statement and presents background on SOS polynomials. 
In Section II VI we define the class of positive operators which we 
utilize. Section 0 provides a controller synthesis condition and 
related inequalities which are later used to prove the main result. 
In Sections IVII we provide the main results wherein we construct 
output feedback controllers. Section I VIII provides the numerical 
results for an example PDE. 

II. NOTATION 

R mxn denotes the set of real m-by-n matrices. S 71 C R nX " 
is the subspace of symmetric matrices. /„ is the identity matrix 
of dimension n x n and we denote I — I n when n is clear 
from context. For any Q C R, C' m (f2) is the space of m-times 
continuously differentiable functions defined on fl. Similarly, for 


AT(a:,£) £ < x 
K 2 (x,£) £ > x 


any fii, Q 2 C R, C ,m,n (ffi x fii) is the space of functions which 
are m-times and n-times continuously differentiable on Hi and 
H 2 respectively. The shorthand u x denotes the partial derivative 
of u with respect to independent variable x. We use L 2 ( 0 , 1) to 
denote the Hilbert space of square integrable functions from [0, 1] 
to R. Unless otherwise indicated, (■,•} denotes the inner product 
on 1 , 2 ( 0 , 1 ) and [| • [| = [| ■ ||x, 2 ( 0 , 1 ) denotes the norm induced 
by the inner product. Similarly, 1 , 2 ( 0 , 00 ; 1 , 2 ( 0 ,1)) denotes the 
Hilbert space of square integrable functions from [0, 00 ) to 1 , 2 ( 0 , 1) 
equipped with the norm 

||/IU 2 (0,oc,;L 2 (0,1)) = ll/(-, i)l| 2 d£^ ■ 

H n ( 0,1) := {y G L 2 : £ L 2 (0,l),i = 1 ,- -n} is 

the Sobolev subspace equipped with inner product (x,y) Hn = 
YlZi=o IF^) ■ F° r Hilbert spaces X and Y, the set C(X , Y) 

is the Banach space of bounded linear operators from X to Y 
endowed with the induced norm || ■ ||£. X denotes the identity 
operator. We define Zd{x) to be the column vector of all monomials 
in variables x of degree d or less. For brevity, we sometimes use 
Z d (x,€) = Z d ([x-,t ]). 

III. PROBLEM STATEMENT 
In this paper, we consider the following scalar parabolic PDE 

wt(x, t ) = a(x)w xx (x, t) + b(x)w x (x, t ) + c(x)w(x, t ) + f(x, t), 

( 1 ) 

where x G [0,1], t > 0, with mixed boundary conditions of the 
form 

«: x ((). t) — 0 . w{l,t)=u(t), (2) 

Here a, b and c are polynomials with a(x) > a > 0, for x G 
[0,1], Additionally, / G L 2 (0, 00 ; £ 2 ( 0 ,1)) is the exogenous input 
and u(t ) is the control input. The output of the system is y(t) = 
w x (l,t). Note that we have also considered several other types of 
observer-controller boundary conditions, which will be listed in the 
section on numerical results. The first goal of the paper is to find a 
control operator T G C (H 2 { 0,1), R) such that if u{t) = Tw{-, t), 
then the closed-loop PDE system is stable. 

Next, using the Luenberger framework, we assume our observer 
has the form 

w t (x,t) =a(x)w xx (x,t ) + b(x)w x (x,t) + c(x)w(x,t) 

+ Oi(x) (y(t) - y(t)), (3) 

with boundary conditions 

ui x (0,t) = 0, w(l,t) = u(t) + 0 2 (y(t)-y(t)) , (4) 

where the function Oi(x) and scalar O 2 must be chosen such that 
the dynamics of the error e(x,t) = w{x, t)—w(x, t ) are stable. The 
second goal of the paper, then, is to find such Oi(x) and O 2 and 
show that if u(t) = t), then the coupled system of parabolic 

PDEs is stable and satisfies 

IMU 2 (0,oo;L 2 (0,1)) < 7II/IU 2 (0,oo;L 2 (0,1))- 

for some 7 > 0. Note that for the system and the observer, 
we assume the existence of classical solutions belonging to 
C 1,2 ((0, 00 ) x [0,1]). This assumption can be validated using the 
analysis presented in [ 2 ] and [ 6 ]. 

A. SOS and Operators 

SOS is an approach to the optimization of positive polynomial 
variables. Given a polynomial f(y), y G R n , the feasibility problem 


of determining if the polynomial is globally positive (f(y) > 0 
for all y G R n ) is NP-hard [3]. To overcome this difficulty, there 
are a number of sufficient conditions for polynomial positivity. A 
particularly important such condition is that the polynomial, p, 
be a Sum-of-Squares (SOS), so that p(x) = X^=i 9 t{x) 2 . for 
polynomials Qi and which is denoted p G E s . The importance of the 
SOS condition lies in the fact that it can be readily enforced using 
LMIs. This is due to the easily proven fact that for a polynomial 
p of degree 2d, p G £ s if and only if p = Z(x) T QZ{x) for some 
Q > 0, where Z(x) is the vector of monomials of degree d or 
less [16]. A recent survey for alternatives to SOS based methods 
may be found in [ 10 ]. 

We can use SOS to construct positive operators on L 2 (0, 1). For 
example, define the operator Vz(x) = M(x)z(x), z G 1 , 2 ( 0 ,1), 
where M is a polynomial. If, for t > 0, M(x) - tGS s , then the 
operator V is positive on 1 , 2 ( 0 ,1). Therefore, we may conclude 
that V is positive on 1 , 2 ( 0 ,1) if there exists a Q > 0 such that 
M(x) — t = Z(x) t QZ(x). By equating the coefficients on the 
left and right-hand sides, we obtain an LMI test for positivity of 
the operator. Of course, the operators considered in this paper are 
significantly more complicated than Vz. 


IV. POSITIVE OPERATORS ON L 2 (0,1) 


In this paper, our results are expressed as optimization over a set 
of positive operators. To solve these optimization problems, we use 
positive matrices to parameterize a subset of positive operators on 
1 , 2 ( 0 ,1) as described in [18]. Specifically, we consider operators 
of the form 

{Vz)(x) = M(x)z(x)+ [ K(x,£)z(£)d£, z G L 2 (0, 1), (5) 
Jo 

with semi-separable kernel 


K(x,0 


Ki{x,0 £,<x 
K 2 (x,£) £ > x’ 


where M : [0,1] —► R and K\, K 2 : [0,1] x [0,1] —^ R are 
polynomials. In [19], we gave necessary and sufficient conditions 
for positivity of multiplier and integral operators of similar form 
using pointwise constraints on the functions M, K\ and K 2 . 
Recently, in [18], these conditions were sharpened - See Theorem 1. 
The following theorem is an extension of this result. 

Theorem 1: Let 


M{x) =Z 1 {x) T U 11 Z 1 {x), 

AT ( x , 0 =Zi {x) T U 1 2 Z 2 (*,£) + Z 2 (£, x) t U 3 1 Zi (£) 

+ [ Z 2 {v,x) T U3 3 Z2{ri,(,)dri 

Jo 

/ x 

Z 2 (p, x) T U 3 2Z 2 {ri, £)dri 
+ f Z 2 {ri,x) T U22Z 2 {ri,^)dp, 

J X 

where K 2 (x,£) = RT(£,®), Zi(x) = Z dl { x) and Z 2 (x, y) = 
Z d2 (x,y) and 


' Un 

c/12 

U 13 ' 


ei/ 

0 

o' 

U21 

c/22 

c/23 

> 

0 

0 

0 

c/31 

c/32 

c/33. 


0 

0 

0 


Then the operator V defined in Eqn. 0 is self-adjoint and satisfies 

£ i|N | 2 < {Vz, z) < e 2 ||a|| 2 , for all z G 1 , 2 ( 0 , 1). 






where £2 = (9 1 + 02)A max ([7), A max ([7) is the maximum eigen¬ 
value of U, and 


9 1 


02 


sup Z\{x) T Z\{x), 
^e[o,i] 


sup 

(*,£)£[0,l]x[0,l] 


/' 


/. 


g(x,€,n)dv+ / g{x,£,v)dri 


g(x,£,v) = Z 2 (r),x) T Z 2 (r),£). 

Proof: The proof is based on the result in [18] and is omitted 
for brevity. ■ 

For convenience, we define the set of multipliers and kernels which 
satisfy Theorem Q] 


H{ dl ,d 2 , ei , e2 } — {M, AT, AT : M,K U K 2 satisfy 
Theorem Q] for di, d 2 , ei, e 2 .} 


Of course, since such operators are positive definite and bounded on 
1/2(0,1), the inverse of these operators exist and are bounded [ 11J. 
However, as will become apparent in subsequent sections, we need a 
method of constructing the inverse of operators defined by elements 
of ~{d 1 ,d 2 ,ei,e 2 }- Fortunately, such methods do exist in literature 
and we use one such method. Using the terminology presented in [8] 
it can be shown that the operators defined by E{d 1 , d2 ,e 1 ,e 2 } are 
the input-output maps of well-posed Linear Time Varying (LTV) 
systems. For this class of operators, the inverse can be constructed 
as explained in [8], 


V. PRELIMNARY INEQUALITIES 
In this section we provide a couple of inequalities which we will 
use for the controller and observer synthesis. We begin by defining 
the operator A : IT 2 (0,1) —>• 1,2(0,1) (infinitesimal genearator) 
which defines the class of PDEs under consideration. 

A=a ^^ +b( ' x ^ +c ^’ (7) 
where recall a, b and c are polynomial functions and a(x) > a > 0, 
for x G [0,1]. Before presenting the inequalities, we define a pair of 
mappings which relate the functions M, K \, K 2 to the derivative of 
the Lyapunov function V = (w,Vw ). The first mapping considers 
(AVz,z) + (z,AVz). 

Definition 1: For scalar £1 > 0 and polynomials a, b 

and c which define the PDE under consideration, we say 

{To, Ti, T 2 , T 3 , T 4 , r 5 , Te} = M ei (AT, AT, AT) if 


T 0 (x) = ( a xx (x ) - b x (x))M(x) + b(x)M x (x ) 


+ a(x)M xx (x) + 2 c(x)M(x) ——aei 


+ a(x) 


2— [ATM-ATM] 


€=x 


Ti(x,£) =a(x)Ki tXX (x,£) + b(x)Ki tX (x,£) 

+ o ( 0^ 1 ,€£(*,0 + 6 ( 0 ^ 1 , e(*.0 

+ (c(x) + c(0) Ki(x,{), 


T 2 (x,£) *=Ti(£, x), 

7T 2 

T 3 = MO) - 6(0)) M{ 0) + a(0)M(0) - —at u 

Ti{ x) =2a(0)K 2lX (0,x) + 7r 2 aei, 

T 5 = (6(1) - 0,(1)) AT( 1) + o(1)M b (1), 

T 6 =2o(l)M(l), 

Ki <x {l,x) — [AI,® (x , £)|x=l]^ =a > ■ 

The second mapping relates the functions N, Pi, P 2 to the deriva¬ 
tive of the dual functional (Aw,Sw) + (SAw,w). 


Definition 2: Given scalar ei > 0 and polynomials a, b 
and c which define the PDE under consideration, we say 
{Qo, Qi, Q 2 , Qi, Q-i, Qs, <26, Qi, Qs} = J\T ei (N, Pi, P 2 ) if the 
following hold 


Qo{x) = 


d_ 

dx 

+ 2 


d_ 

dx 


a(x)N(x) — b(x)N(x) 


-7^ l a ( x ) (A(X,C) - Plfa: 0)1 


J £=x 


+ 2 N(x)c(x) ——aei, 


Q i(M) = 


dx 


dx 


d 

+ d£ 


a{x)Pi(x,£,) - b(x)Pi(x,£) 

^a(0Pi(x,0-K0Pi(x^) 


+ (c(x) + c(Z)) Pi{x,Z), 

Q 2 (at,0 =Qi{f,x), 

7T 2 

Q 3 = (0,(0) - 6(0)) N{ 0) + o(0)AT(0) - —aei, 

Qi(x) =2 (a*(0) - 6(0)) P2(0,®) + 2a(0)P 2 , x (0, x) 

2 

+ 7r aei , 

Q 5 =(6(1) - a x (l))N(l) - a(l)AT(l), 

Qe(x) =2(6(1) - a a (l))Pi(l, x) - 2a(l)Pi,*(l, x), 

Qr =2a(l)N(l), 

Qs(x) =2a(l)Pi(l, x), 

Pl,x(l, x) = [Pl,ar(®,0|s=l]£=* ■ 

The proofs of the following lemmas are provided in the appendix. 
The first allows us to represent V = {AVz, z) + {z, AVz) 


Lemma 1: For any {AT, AT, AT} € E { d 1 ,d 2 ,e 2 ,e 2 }, 0 < ei < 
£2 < 00, let {To, 7T, T 2 , T 3 , TT, T5, Tq} = M tl (M,Ki,K 2 ). 
Then, for any w G T2(0,1), if the operator V is given by 

(Vw) (x) = M{x)w(x) + ( K(x,£)w(t;)dx, (8) 

Jo 

with 

K(x,0 = 

and operator A is given by Equation 0. we have that 


Ki{x,i) £ < x 
K 2 (x,f) £ > x 


{AVz, z) + (2, AVz) 

< (Tz, z) + 2(0) ^T 3 z(0) + J Td{x)z{x)dx S j 
+ z(l)(T S z(l)+T 6 Z x (l)), 

where 2 = V~ 1 w for any w G H 2 { 0,1) with w x (0) = 0. Here we 
define the operator T as 

(Tz) (x) = T 0 (x)z(x) + [ T(x,€)z(£)d£, 

Jo 

with 

T(x,0 = 

The second lemma allows us to represent the derivative of 

(Aw, Sw) + (SAw, w). 

Lemma 2: For any {N, Pi,P 2 } G ^{d 1 ,d 2 ,e 2l e 2 }, 0 < 

ei < £2 < 00, let {Qo,Qi, Q 2 ,Q 3 ,Qa,Q 5 ,Qo,Qi, Qs} = 
Af ei (N, Pi, P 2 ). Then, for any 2 G L 2 (0, 1), if operator S is given 


Ti(x,fi) i<x 
T 2 (x,T) £>x' 














by 


1 


(Sz)(x) = N(x)z(x) + P(x,t,)z(£)dx, 


(9) 


with 


and operator *4 is given by Equation 0, we have that 

(>4u>, <S«t) + w) 

< (Qw,w) + w(0) ^Q3 w( 0) + j Q±(x)w(x)dx 
+ w(l) ^Qsw(l) + J Qo(x)w(x)dx^ 

+ 10^(1) ^Qrw(l) + J Qs(x)w(x)dx' S j , 

for any w £ H 2 ( 0,1) with w x ( 0) = 0. Here we define the operator 
Q , for any z £ £2(0, 1), as 


(Qz) (x) = Q 0 (x)z{x) + / Q(a;,f)z(f)d£, 


with 


Q(*,£) = 


k Q2{x,0 £>x 

VI. OUTPUT FEEDBACK CONTROLLER SYNTHESIS 


Our approach to design of output-feedback controllers is based 
on three steps. First, we design the control operator T which maps 
that state to the control input as u(t) = J-w. However, because we 
cannot measure the state, we find function Oi(*) and scalar O 2 
which define the observer which outputs an estimate of the state 
w. Finally, we prove that the controller coupled to the observer as 
u(t) = J-w produces a closed-loop system with bounded L 2 gain 
from exogenous input to controlled output. 


A. Control Design 

We begin by designing the control operator T £ 
C. (H 2 ( 0,1),R). Consider the following observer dynamics 


w t (x,t) =a(x)w xx (x,t) + b(x)w x (x,t) + c{x)w{x,t) 

+ Oi(x)e x (l,t), (10) 

with boundary conditions 

wv(0,f) = 0, w(l, t) = u(t) + C> 2 e x (l,t). (11) 

The following lemma defines the operator T . 

Lemma 3: Suppose there exist {M, Ki, K 2 } £ Hd 1 ,d 2 ,ei,£ 2 an d 
Ti such that {To, Ti, T 2 , T 3 , T 4 , T 5 , T&\ — M E1 (M, Ki, K 2 ) and 

T3 < 0, T 4 ,( x ) = 0. 

Let u{t ) = where T — ZV~ 1 , V is as in Eqn. 0 and 

the operator Z is defined as 

(■ Zg ) (x) := Z\g x {l) + f K 1 {l,x)g(x)dx. 

Jo 

where Z\ is any scalar such that T 5 Z 1 > —TqM{ 1). Now, if 
V(w(-,t)) = {w(-, t), t)) where V is as in Equation 0. 

Then for any solution w of Eqns. m and OB with input e x , we 
have that 

^V(w(-,t)) < (Tz(-,t),z(-,t)) + 2{0 1 (-)e x {l,t),z(-,t)) 


— gz(l, t) 2 - ?-^-z(l,t)e x (l,t), 

for some g > 0 where z = V~ 1 w and T is defined in Lemma [I] 
Proof: We begin by taking the time derivative of the Lyapunov 
function V 0 (w(-,t)) along the trajectories of fllO Ml 11 1 

= (Aw(-,t),'P~ 1 w{-,t)) + (P" 1 w(-,f),^lr 6 (-,f)} 

+ 2 ( Oi(-,t)e x (l,t),'P~ 1 w(-,t )) , 

where we use that V^ 1 is self-adjoint have simplified the derivative 
using the definition of operator A provided in Equation 0. We 
rewrite this as 

= (AVV^w{-, f), P _ 1 u>(-, t)) + (P _ 1 w(-, t), AVV^wi-, t)) 
+ 2 (Oi(-, t)e x (l, t)) . 

Now define z = V~ 1 w, then 

^I4(w(-,f)) =(AVz{-,t),z(-,t)) + (z{-,t),AVz(-,t)} 

+ 2 (Oi(-, t)e x (l, t)). 

Now, applying Lemma Q] and using the facts that T 3 < 0 and 
Ti{x) = 0 produces 

~ffV 0 {w{-, t)) < (T z(-, t ), £(■, t )) + 2 (Oi(-, t)e x { 1 , t), z(-, t)> 

+ z(l, f) (T 5 z(1, t) + T 6 z x ( 1, t )). (12) 

Since z = V~ 1 u>, w = Vz. Thus 

= M(l)z(l,t) + [ Ki(l,x)z(x,t)dx. (13) 
Jo 

From the boundary condition in 03 we get 

w(l, t) = u(t) + C> 2 e x (l,t) =J r w(-,t) + 02 e x (l,f) 

=J r VV~ 1 Ul(-,t) + 02ea:(l,t) 
=Zz(-,t) + 0 2 e x (l, t). 


Using the definition of Z, 

w(l,t) = Ziz x (l,t) + f Z 2 {x)z(x,t)dx + C> 2 e x (l,t)- 
Jo 

Substituting into Equation 03 and using the definition Z 2 (x ) = 
AT(l,a:) we get 

- ^re x (l,t)- 

Substituting this expression into 03 

d V 0 (w(-, t)) < (Tz(-, t), £(■, t)) + 2 (Oi(-, t)e x {l, t), £(■, t)) 


dt 


+ (J5 + 

T0O2 


TeM( 1) \ 


Zi 


Z 1 J 

z(l,t)e x (l,t)- 


(1 A? 


Now, since Z\ < 0 is a scalar such that T5 + T§M(X)/Z\ < 0, 
there exists a scalar g > 0 such that T5 + T§M(1)/Z\ = —g. 
Hence 
d . 


dt 


Vo (w(-,t)) < (Tz(-, t), z(-, t)) +2{Oi{-)e x {l,t),z(-,t)) 


— gz(l, t) 2 - 


-z(l,t)e x (l,t)- 


Z 1 







B. Observer Design 

We now design the function 0\{x) and scalar O 2 which define 
the observer. We begin by subtracting Equations ©-© from 0- 
0 to obtain the dynamics of the error variable e = w — w given 
by 

e t {x,t) =a(x)e xx (x,t ) + b(x)e x (x,t) + c(x)e(x,t) 

+ Oi{x)e x (l,t) - f{x,t), (14) 

with boundary conditions 

e x (0, f) — 0, e(l, f) = C^e^l, f), (15) 

where we have used the definition of the measurement y(t) = 
w x (l,t) and y(t) = w x (l,t). We present the following lemma. 

Lemma 4: Suppose there exist {N, Pi, P 2 } G ’ELd 1 ,d 2 ,ei,e 2 ’ suc h 
that 

Q3 < 0 , Qi{x) = 0 , 

where {Qo, Qi , Q2, Q3, Qi, Qs, Q 6 , Qi, Qs} = J\Te 1 (N, Pi, P2) 
(See Defn.©. Let S be defined as in Eqn. 0. Then choose a scalar 
O2 < 0 such that 

Qs + 77-Q7 < 0, 

U 2 

and let Oi(x) = {S~ 1 Ri) ( x ) where Ri(x) = 

— \ (C> 2 Q 6 (x) + Qs{x)). Define V e (e) = (e,Se). Then for 

any e, f which satisfies Eqns. cs-cnDi with Oi(x) and O 2 as 
defined here, we have 

^V' e (e(-,f)) <{Qe{-,t),e(-,t)) + 2 (f(-,t),Se(-,t)) 
-Ce(l,f) 2 , 

for some scalar f > 0 where Q is as defined in Lemma [2] 

Proof: We begin by taking the time derivative of the Lyapunov 
function V e (e(-,t)) along the trajectories of <14M15t . yielding 

^14(e(-,f)) = (e t (-,t),Se(-,t)) + (e(-,t),Se t {-,t)) 

= {Ae(-,t)),Se(-,t)) + (e(-,t)),SAe(-,t)} 

+ 2 (e(-, t), (SOi) (-)e x (l,t)) + 2 {f(-,t),Se(-,t)) , (16) 

where we have again used the definition of A from Eqn. 0 and 
we have also usd the fact that S is self-adjoint. Now, since from the 
theorem statement we have that Q 3 < 0 and Q 4(2) = 0, applying 
Lemma [2] produces 

~fjV e { e {-,t)) < (Qe(-, t),e(-,t)) + 2 (/(•, t),Se(-, t)) 

+ e(l,t) ^Qse(l,t) + j Q&{x)e(x,t)dx^j 

+ e x (l,t) (^Q 7 e(l,t) + I! Qs{x)e(x,t)dx' S j 

+ 2(e(-,t),Ri(-)e x (l,t)), (17) 

where we have used the fact that since Oi(o;) = (.S^-Ri) (*), 
Ri{x) = (SOi) (x). We have the boundary condition e(l,t) = 
02 e x (l,t) and since O 2 < 0, we have that e x (l,t) = e(l,f)/02- 
Substituting in fT7t . 

< (Qe(-,t),e(-,t)) +2 (f(-,t),Se{-,t)) 

+ (qs + e ( i >^) 2 


+ e(li t) 


n 


1 2 \ 

Qg{x) + -pr-Q&{x) + -f-Ri(x) e(x,t)dx. 
U 2 U 2 / 

( 18 ) 


Since O 2 < 0 is a scalar such that Qs + Q7/O2 < 0, let 

C = —(Qs + 09) 

Then £ > 0. Now, using the definition of Ri(x) we get that 

Qe{x) +-^-Qa(x) +-^-Ri(x) = 0. (20) 

u 2 U 2 

Substituting Eqns. (1 1 9l>-ll20b into Eqn. ( 1181 . we find 


d 

dt 


K(e(-,f)) <(Qe(-,t),e(-,t)) + 2 (f(-,t),Se(-,t)) 

-Ce(l ,tf. 


C. Output Feedback Based Control 

We now have the following set of coupled parabolic PDEs. 

w t (x,t) =a(x)w xx (x,t) + b{x)w x {x,t) + c(x)w{x,t) + f(x,t), 

( 21 ) 

Wt(x,t) =a(x)w xx (x,t) + b(x)w x (x,t) + c(x)w(x,t) 

+ Oi(x) (w x (l,t) - w x (l,t)), (22) 

with boundary conditions 

^7,(0, i) = 0 , w(l,t) = Tw{-,t), ( 23 ) 

thx(0, f) = 0, ti)(l,f) = Fw(-,t) + O2 (w x (l,t) - w x (l,t)). 

( 24 ) 

We now prove that the previously designed controller and the 
observer can be coupled such that norm of the system state remains 
bounded in the presence of an exogenous input. 

Theorem 2 : Suppose there exist scalars 0 < ei < £2 < 00, 
5 , /3 > 0 d\, d2 G N and polynomials {N, Pi,P2} G Ed li d 2ieijE2 
and {N, Pi, P 2 } G Ed 1 ,d 2 , ei ,t 2 , such that 

{—To — 2 5 M, —Ti — 2 SKi, —T2 — 2<5A'2} G Ed 1 ,d 2 ,o,p, 
{—Qo — 25 N, —Qi — 25 Pi, —Q2 — 25P2} G ’S.d 1 ,d 2 ,o,p, 

T 3 < 0, Ti(x) = 0, Qi < 0, Q4(s) = 0, 

where {To, Ti, T2, T3, T4, T5, To} = A 4 ei (M, Ki, K2) and 
{Qo, Q 1, Q21 Q35 Qi: Q 5: Qg: Qt ; Qs} — ^e\ (-(V5 Ri 5 R2 ) as pro¬ 
vided in Definitions Q] and [ 2 ] respectively. Let V be defined as in 
Equation 0 and S as in Equation 0. 

Then for any solution \w{x, t ) w{x, f)] of the coupled dynam¬ 
ics { 2 D-{ 24 ]i. if J- is given by Lemma [ 3 ] and Oi {x) and O2 are 
given by Lemma [ 4 ] there exists a scalar 7 > 0 such that 

|M|l 2 (0,oo;L 2 (0,1)) — 7ll/l|i 2 (0,oo;I, 2 (0,l))i 

for any / G I/2 ( 0 , 00; 1/2(0, 1 )). 

Proof: For the Lyapunov function V 0 (w(-,t)) = 

we have from Lemma [ 3 ] that there exists 

scalar p > 0 such that 

^y o {w{-,f)) < ( Tz(-,t),z(-,t )) + 2 (Oi(-)e x (l, t),z{-,f)) 

- pz(l,t) 2 - ^-^-z{l,t)e x {l,t). 

We have from Lemma | 4 ] that e x (l,t) = e(l,t)/02. Therefore 
^-Vo{w{-,f)) < (Tz(-,t),z(-,t)) + J- <Oi(-)e(l,f),z(-,f)} 



(32) 




(25) 


For the Lyapunov function I4(e(-,i)) = (e(-,t),Se(-,t)), we 
have from Lemma [4] that there exists a scalar £ > 0 such that 

^K(e(-,i)) <{Qe(-,t),e(-,t)) + 2(f(-,t),Se(-,t)) 

-Ce(l ,t) 2 . 

From Equations ( I25b - d26b we conclude that for any A > 0 

d ,, , „, .. , d 


(26) 


d V 0 (u,(,t)) + A-VM,t)) 

< A(Qe(-,t),e{-,t))+ 2A (f(-,t),Se(-,t)) 

+ < 


’«(•,*)" 


T 0 

o ' 



z{l,t) 

) 

-k —fl 

t 6 

2Z, 


z(l,t) 

_e(l,f)_ 


k k 

\ 


_e(l,f)_ 


(27) 


where ( Oy ) (x) = ^Oi(x)y(x), for any y G 1/2(0,1), and the 
inner product is defined on 1/2(0, 1) x 1/2(0,1) x 1/2(0,1). Now, 
since {—To — 28M, —T\ — 28K\, — T 2 — 25 A 2 } G 'Sd 1 ,d 2 ,o,p, we 
have that T + 28V < 0. Therefore, for any 0 < 9 < 8 , it can be 
established using Schur complement that for a large enough A > 0, 


T + 2 ev 

~k 

~k 


0 


o 

Ta 

~ “ 2 Z-l 

* -A< 


< 0 . 


Therefore 


T 0 


O 


★ —u — 26s_ 

P 2Z\ 

~k ★ —A(( 


< 


-2 ev 0 0 

* 00 

★ ★ 0 


Substituting into Equation J27b . we get 

±V 0 W,t)) + Af t VM,t)) 

< A (Qe(-,t),e(-,t)) + 2A 
- 26 (z(-,t),Vz(-,t)). 

Let V 0 (w(-,t)) + AV e (e(-,t)) = V(t), thus 
d . 


dt 


V ( t ) + 20 (z(-, t),Vz(-, t)) <A (Qe(-, f), e(-, t)) 

+ 2A </(-,t),«Se(-,i)). 


Adding A5 (5e(-, t)), e(-, f)) — y (/(•, t), Sf(-, t)} to both sides, 
jV{t) + 26 (z(;t),Vz(; t)> + Ai (<Se(-,f)), e(-, t)> (28) 


--</(-, t),5/(-,t)) 


< A 


e (', 0 



5 


Q + <55 

Vr 


-fo L . 


e(-,f)' 

/M) 


(29) 

(30) 


Since {—Qo — 28N,—Qi — 2<5Pi,— Q 2 — 28 P 2 } G Sd ll( j 2l o,/3, 
we have that Q + 2<5<S < 0. Flence, using Schur complement we 
conclude that 

'Q + 8 S S' 

* -is 


< 0 . 


Therefore, from Equation 130b we conclude that 

^ V'(i) + 26 (z(-, t),Vz(-, t)) + A 8 (,Se (•, t)), e(-, f)) 
<j{f(;t),Sf(-,t)). (31) 

Since the operator S is defined using {N, Pi, P 2 } G 


we have from Theorem |T| that, for all g G T2(0,1), 
eillsll 2 < ( Sg,g) < t 2 \\g\\ 2 . 

Similarly, since V is defined using {M, Ki,K 2 } G ^d 1 ,d 2 ,e 1 ,e 2 , 
using Theorem [T] it can established that 

-||S|| 2 <(P- 1 5,3>< 1 II<?I| 2 - 

£2 £l 

Since (z(-,t),Vz(-,t)) = (V~ 1 w(-,t),w(-,t)'), from the previous 
expression we have that 


_l_ 

£2 


\w{-,t )\\ 2 < (V = {z(-,t),Vz(-,t)). (33) 


Substituting Equation ( 133b in Equation 13 1 b and using 132b . we get 

±V(t) + 2^\\w(.,t )\\ 2 + A 8 ei\\e(;t )\\ 2 < ^||/(-,f)|| 2 . 
at 62 0 

Integrating in time from 0 to some 0 < T < oo, we get 

V(T)-V(0) + 2— [ T \\w(;t)\\ 2 dt + A 8 ei f T \\e(;t)\\ 2 dt 

£ 2 Jo Jo 

< ^ l^\\f(;t)\\ 2 dt. (34) 

Now, V(T) = V 0 (w(-,T)) + AV e (e(-, T)) > 0. Additionally, if 
we assume zero initial conditions, then U(0) = I4(ui(-,0)) + 
AV e (e(-, 0)) = 0. Therefore we conclude from Equation 134b that 


f 

Jo 


\\w(;t)\\ 2 dt+ ||e(- 


/ 


ll’ d 's 1 ? 


f 

Jo 


ll/(’i t)\\dt, 


where v = min |2^-, A<5ei|. Since / G L 2 (Q, oo; T2(0,1)), 
taking the limit T —>• oo, we get 


IHli 2 (0,oo;i 2 (0,l)) + ll e l| 2 2 (0,oo;I, 2 (0,l)) 
- -^7ll/lll 2 (0,oo;L 2 (0,l))- 
Hence, we conclude that 


IH|,L 2 (0,oo;£ 2 (0,1)) <A 
||e||i 2 ( 0 ,oo;i 2 ( 0 ,l)) <]J-J^-\ 


l|i 2 (0,oo;L 2 (0,l)) 5 

|U 2 (0,oo;L 2 (0,1))- 


Since e = w — w, w = w — e. Therefore 

ll w |U 2 (0,oo;Z, 2 (0,l)) < ||w’||z, 2 (0,oo;L 2 (0,1)) + II e IU 2 (0,oo;Z, 2 (0,1)) 

“ ^^ 11/11^(0,00^2(0,1))- 

Setting 7 = 2 completes the proof. ■ 

VII. NUMERICAL RESULTS 

In this section we consider a couple of examples on which we test 
the conditions of Theorem [2] using SOS and SDP. These numerical 
results are obtained using the Matlab toolbox SOSTOOLS [20]. 

We consider the following two PDEs. First consider the classical 
heat equation with an unsteady source term. 

w t (x,t) = w xx (x,t) + Xw(x, t) + f(x,t), (35) 

Without feedback, this system is unstable for A > iv 2 /4. Next, we 
consider a randomly generated PDE. 

w t (x,t) = (x 3 — a: 2 + 2) w xx (x,t)+ (3* 2 - 2x) w x (x, t) 

+ (—0.5® 3 + 1.3i 2 — 1.5x + 0.7a: + A) w(x, t) + f(x, t ) 

(36) 






















TABLE I: Maximum A as a function of polynomial degree di = 
d ,2 = d for which we can construct output feedback boundary 
controllers for PDE ( f35l) . 


d= 7 8 

9 

10 11 

A = 12.69 16.01 

17.96 

17.96 21.97 

TABLE II: Maximum A as a 
d 2 = d for which we can 
controllers for PDE 136b. 

function of polynomial degree di = 
construct output feedback boundary 

II 

"C3 

6 

7 8 

A = 18.75 28.78 

32.03 

32.03 39.16 


with boundary conditions 

ua(M) = 0, w(l ,t) = u(t), (37) 

By using stability analysis and numerical simulation, we estimate 
that PDE is unstable for A > 4.66. 

In these examples, we find the maximum A, using a bisection 
search, for which we can construct stabilizing output-based bound¬ 
ary feedback controllers. We test the conditions of Theorem [2] with 
ei = 0 . 001 , £2 = 1, <5 = 0.001 and increasing values of d\ 
and d 2 - Table [T] presents the maximum A > 0 for which we can 
construct output feedback controllers for PDE d35b as a function 
of the degree di = d 2 = d of the polynomials which define the 
controller, observer, and Lyapunov function. Table |TT] presents the 
maximum A > 0 for PDE 136b . 

The numerical results suggest that increasing the degree di = 
d 2 = d of the polynomial representation leads to higher values of 
A > 0. Moreover, the value of A > 0 does not appear to be upper 
bounded, which implies that the method is asymptotically accurate. 
That is given any A > 0, we conjecture that we can construct output 
feedback controllers for a large enough degree d. 

Figures 11121 represent the simulation of PDE ( 136b with A = 39 
subject to the output feedback based control in the presence of 
exogenous input f(x, t ) = e -i cos(7rf) (1 + sin(0.l7ra;)). 


o.i 



Fig. 1: State of PDE ( 136b with point observation and point actuation. 

One of the key technical advances of this paper is the use of 
semi-separable kernels K i, K 2 , Pi and P 2 and this advance leads 
to significantly more complex stability conditions. Therefore we 
wish to establish if the inclusion of the variables K 1 , K 2 , Pi 
and P 2 does, in fact, provide any significant performance gain. 
In order to do this, we check the conditions of Theorem [2] while 
setting A'i = K 2 = Pi = P 2 = 0 (similar to our previous 
approach [7]) and applied these conditions to the example PDEs. 
Table Hill presents these results for PDE ( 135b and Table HVl presents 
results for PDE C36J- 



Fig. 2: Boundary control effort w(l, t ) = u(f) for PDE ( f36t . 

TABLE III: Maximum A as a function of polynomial degree, di = 
d 2 = d for PDE ( 135b for which we can construct controllers using 
with Ki = K2 = Pi = P2 = 0. 

d = 1 2 3 4... 10 

A = 3.91 4378 4~88 4~88 


Comparing Tables imnvl with Tables HE] we observe that the 
inclusion of kernels K 1 , K 2 , Pi and P 2 allows the construction of 
output feedback based controllers for significantly higher values of 
A. Moreover, by setting I\ 1 = I\ 2 = P 2 = P 2 = 0, the numerical 
results appear to show an upper bound to the A for which we can 
design controllers without the use of these kernel functions. We 
conjecture, therefore, that kernel functions are a necessary part of 
any Lyapunov-based method for analysis and control of PDEs. 

Finally, as we previously stated, the SOS conditions for the 
design of output feedback controllers can be easily modified for 
systems with other types of boundary conditions. To this end, we 
provide the numerical results for controller synthesis for PDEs ( 135b 
and ( 136b with boundary conditions defined in Table [V] Table 1 VII 
illustrates the maximum A for which we can construct output 
feedback controllers as a function of di = d 2 = d for PDE ( 135b 
with boundary conditions and outputs given in Table m Similarly, 
Table IVIll illustrates these results for PDE 136b . 

We note that the backstepping method has been applied to 
Example 135b and is also able to construct exponentially stabilizing 
output feedback boundary controllers for arbitrary A > 0 (see [12]). 
Therefore, while we cannot necessarily claim any improvement in 

TABLE IV: Maximum A as a function of polynomial degree, di = 
d 2 = d for PDE 136b for which we can construct controllers using 


with Ki = K2 = 

Pi = P2 = 0. 


d 

= 1 2 3 4... 

10 

A 

= 3.51 5.47 6.64 6.64 

TABLE V: Alternative boundary conditions 
PDE 0D-<36). 

and outputs 


Boundary Condition 

Output y{t) 

Dirichlet 

tn(0, f) = 0 
w(l, t ) = u(t ) 

W X { 1 , t ) 

Neumann 

w x (0,t) = 0 
w x (l,t) = u(t) 

w(l,t) 

Robin 

w(0, t) + Wx (0, f) = 0 
w(l,t) + = u(t) 




















TABLE VI: Maximum A as a function of polynomial degree, d\ = where 

d ,2 = d for PDE ( 135b with boundary conditions and outputs given 

in Table 0 L, = 



d= 7 

8 

9 

10 

11 

r 2 

Dirichlet 

A = 14.25 

17.96 

17.96 

24.21 

25.78 

Neumann 

12.69 

16.01 

17.96 

17.96 

21.97 

r 3 

Robin 

11.71 

14.45 

16.40 

17.96 

18.84 

TABLE VII: Maximum A as 

a function of polynomial degree, di = 


d 2 = d for PDE d36b with boundary conditions and outputs given 

r 4 

in Table [V] 







d = 4 

5 

6 

7 

8 


Dirichlet 

A = 21.87 

33.59 

36.71 

36.71 

44.53 


Neumann 

18.75 

29.78 

32.03 

32.03 

39.16 

r 5 

Robin 

14.16 

26.66 

28.90 

28.90 

30.46 



performance over this established methods, our approach is at least 
competitive and may have certain advantages such as relative ease 
of implementation (changing the system is a one-line edit) and the 
fact that our approach does not require numerical integration of a 
PDE. 


f 1 d 2 

J z{x)a{x)-^^[M{x)z[x)]dx, 

■ [ z(x)b(x)-^~ [M{x)z(x)\ dx, 

Jo UX 

lo (/ K i( x ’0 z (0<%) dx 

+ J z ( x ) a ( x )-^ (^J K 2 (x,i)z(i)d^jdx 

L Z ^ b ^'§x (/ dx 

+ lo K ^ x ’^ z ^ d ^j dx ’ 

- f z(x) 2 M(x)c(x)dx 

Jo 

rl nx 

+ / z(x)c(x)Ki(x,£)z(£)dt;dx 

Jo Jo 

+ / z(x)c(x)K 2 (x,f)z(£,)d^dx, 

Jo Jx 

where we have used the fact that 

Ki(x,£) £<x 


K(x,0 = 


K 2 (x,£) t; > X 


VIII. CONCLUSIONS 

In this paper we developed an algorithmic approach for designing 
output feedback boundary controllers for a class of linear scalar 
valued inhomogeneous parabolic PDEs. Our approach is based on 
a parameterization of positive multiplier and integral operators with 
semi-separable kernels. We tested the approach on homogeneous 
and inhomogeneous PDEs using several different types of bound¬ 
ary feedback and several different types of point measurements. 
Furthermore, we tested our approach with and without kernel 
functions to determine if kernel functions are a necessary part of 
Lyapunov theory for PDEs. Our numerical results indicate that 
kernel functions are a necessary part of Lyapunov functions for 
PDEs. Further, our numerical results indicate there is little or no 
conservativity in the method and that our approach is competitive 
with well-established approaches such as backstepping. Note that 
as yet, the observer-based controllers in this paper are not optimal 
in any norm. Therefore, an obvious future direction of this work is 
to extend our approach to Woo-optimal control. 


APPENDIX 


To prove Lemmas Q] and [2] we use the following identity. 
Lemma 5 ([9],[13]): let w G H 2 ( 0,1) be a scalar function. 
Then 


/ 


(w(x) — w(0)) 2 dx < —? [ w x (x) 2 dx. 

7T 2 


Proof: [Lemma U) We begin by considering the following 
decomposition 

{AVz, z) + {z , AVz) 

= 2 lo [(7 ? 2)(k)] + 6 (a:)^ [(' Pz)(x)]J z{x)dx 

+ f c(x){Vz){x)z{x)dx 
Jo 

= 2(ri + L2 + r3 + r4 + r 5 ), (38) 


Before we proceed, we calculate the boundary condition at x = 
0. Since 0 = V~ 1 w, for any w G H 2 ( 0,1) with w x (0 ) = 0, we 
have that w = Vz. Using the definition of V, 

w x (0) = M*(0)*(0) + M(0)z x {0) + [ A' 2 , 40 , x)z{x)dx, 

Jo 

where we have used the fact that K\(x,ff) = K 2 (f,x). Since 
Wa,( 0 ) = 0 , we conclude that 

— M(0)z x {0) = M x (0)z(0) + f K 2 , x {0,x)z(x)dx. (39) 

Jo 

Applying integration by parts twice and using the boundary condi¬ 
tion at x = 0 , we get 

ri=— f z x (x) 2 a(x)M(x)dx 

Jo 

1 f 1 

+ - J z(x) 2 [a xx {x)M{x) + a{x)M xx (x)\ dx 

+ i*(l ) 2 [a(l)Mx(l) - a*(l)Af(l)] 

+ —z(0) [a(0)M x (0) + ax(0)A/(0)] 

+ z(l)a(l)M(l)z x (l) + z(0) [ a(0)K 2 , x {0,x)z(x)dx. 

Jo 


From Theorem Q] it is readily established that M(x) > ei. 
Additionally, we have that a(x) > a. Therefore, a(x)M(x ) > aei 
and we may apply Lemma [5] to produce 

,2 ri 

\ — z( 0)) 2 dx. 


f 


— I z x (x) 2 a(x)M(x)dx < -^-ae 1 


[ (z(*) 

Jo 


Therefore, 
Ti < 


/ ’ z{Tf (I 


9 [a xx (x)M(x) + a(x)M xx (x)\ —— aei ) dx 


+ -z(l) 2 [a(l)M x { 1) - a x (l)M(l)] 


(38) 






+ 2 *( 0) 2 o( 0)M X (0) + a x (0)M(0) ——crei 

+ z( 0 ) J ^a(0)K 2 ,x(0, x) + ~ 2 ae ^j z(x)dx 
+ z(l)a(l)M(l)z x (l). 

Similarly, applying integration by parts once, 

r i 

\ 2 | 


(40) 


r2 4 


z(x) [b(x)M x (x ) — b x (x)M(x )] dx 


o 


+ \z{lfb{l)M(l) - i.~(0) 2 6(0)M(0). 


(41) 


Now, note that for (M, Ki, K 2 ) G Hd l!t i 2iei , e2 , we have 
Ki(x,f) = I\ 2 (^,x) and thus I\i(x,x) = K 2 (x,x). Utilizing 
this property and applying integration by parts twice 


r 1 

f r 

Jo Z ^ 

a{x) 


d_ 

dx 


[K 1 (x, 0 -K 2 {.x,^)} 


J £=X 


dx 


n 1 nx 

+ / z(x)a(x)K 1:XX (x,£)z(£)d£dx 

Jo Jo 

+ / / z(x)a(x)K 2 ,xx(x,£)z(€)d£dx 

Jo Jx 


Dividing the double integrals 

f 1 


r 3 = 


d_ 

dx 


[K 1 (x,f)-K 2 (x, 0] 


[ 

Jo 

+ \J J z(x)a(x)Ki, xx (x,£)z(£)dtda 

+\ii: 


J £=x 


dx 


z(x)a(x)K 2 , xx {x, £)z(£)d£dx 


' 0 J x 

n 1 nx 

Jo Jo 


z(x)a(x)K ltXX (x,f)z(f)d^dx 


+ \ [ [ z(x)a(x)K 2 ,XX {x,0z((.)d£dx 
Z Jo Jx 

Changing the order of integration, switching between x and £ and 
using the fact that = Kz(£,x) in the last two double 

integral produces 


r 3 


r 1 

/ r 

L <x) 

a(x) 


— [K 1 (x^)~K 2 (x,0] 


J £=x 


dx 


Finally, applying a change of order of integration as applied to T 3 
and F 4 , 

= f z(x) 2 M(x)c(x)dx 
Jo 

+ \ j o J o z{x)[c(x)+c(£)]K 1 (x,£)z{£)d£dx 

+ \f J z(x)[c(x) + c(£)\K 2 (x,£)z(£)d£dx (44) 

Substituting d40b -( l44b in d38b and using Definition |T] 

{ AVz,z) + ( z,APz ) 

nl nl nx 

< / z(x) 2 To(x)dx + / / z(x)T\(x,( > )z(f)df y dx 

Jo Jo Jo 

+ [ [ z(x)T 2 (x,(,)z(£,)d£dx 
Jo Jx 

+ z(0) ^T 3 «(0) + J T4(x)z(x)dx^ 

+ z( 1) (T 5 z(l) + T 6 z x (l)). 

Finally, using the definition of operator T, 

(AVz, z) + (z, AVz) 

< {Tz, z) + z(0) ^T 3 z( 0 ) + J T 4 ,{x)z(x)di 
+ 2 ( 1 ) (T 5 z(l) +T 6Zx (l)). 


Proof: [Lemma [2) Using the self-adjointedness of operator S 
we begin with the following decomposition 

(Aw, Sw) + (SAw, w) 

= 2 ( Aw,Sw) 


>/’ 


= 2 / (a(x)w xx (x) + b(x)w x (x) + c(x)w(x)) (Sw)(x)dx 


= 2(r 1 + r 2 + r 3 + r 4 + r 5 ), 

where 


(45) 


+ 


/ z(x) [a(x)K ltXX (x,f) + a(f)K lti( (x,f)]z(^)d^dx 
* Jo Jo 

+ \ J J z(x) [a(x)I\ 2 ,xx(x, £) + a(^)K 2 ,^(x, 0 ] z(£)d£dx. 

(42) 

Applying integration by parts once and employing the previously 
performed change of order of integration 


Ti = / w xx (x)a(x)N(x)w(x)dx, 

Jo 

r 2 = / w x (x)b(x)N(x)w(x)dx, 

Jo 

n 1 nx 

r 3 = / w xx (x)a(x) / Pi(x,£)w(£)d£dx 
Jo Jo 

+ / w xx (x)a(x) / P 2 (x,£)w(£)d£dx, 

Jo Jx 

n 1 nx 

r 4 = / w x (x)b(x) / Pi(x,£)w(£)dfdx 
Jo Jo 

+ / w x (x)b(x) / P 2 (x,£)w(£)d£dx, 
Jo J x 


r 4 


^ Jo Jo z(x)[b(x)K ltX (x,f)+ b(^)K 1 ^(x,^)]z(f)z(x)dfdx 

+ 7 ; J o J z(x) \b(x)I<2,x (a;, £) + b(£)K 2ti (x, 0] z(g)z(x)d£dx. 


( 43 ) 


f 


Ts = I w(x) N(x)c(x)dx 

nl nx 


a: 

a: 


w(x)c(x)P\(x, £)w(£)dt;dx 


w(x)c(x)P 2 (x, £)w(£)dt(dx. 










Here we have used the fact that 

r Pi(x,t) £<x 


P{x,0 = 


Pi{x,i) t>x 


V -/ 1 


Ti = - w x (x)a(x)N(x)dx 


+ 


\j Q ]^l a ( x ) N ( x )] w ( x ) 2dx 


IT < 


2 (^^2 [o(®)-^(*)] - Y aei ) dX 


r z = ~\f w{x) 2 -^ [b(x)N(x)]dx 

+ ~ 6 (l)JV(l)«,(l) a - ±b(0)N(0)w(0) 2 . 


r 3 = / w(x) 


[a{x){Pi(x,£) - P 2 (z,f))] 


(* 


e,=x 


dx 


+ J w{x) J (^^3 [a(a:)Pi(a:,£)]J w(£)d£dx 

+ Jo f (JJ? w(£)d£dx 

-w(l) [ (a x (l)Pi(l,x) + a(l)Pi, x (l,x))w(x)dx 
Jo 

+ w( 0) f (a x (0)P2(0, x) + a( 0 )P 2 ,x( 0 , x)) w(x)dx 
Jo 

+ w x (l) / a(l)Pi(l, #)u;(x)dx. 

Jo 

We can divide the two double integrals as 

■ [ * x) '{[i 


r 3 = 


[o(a;)(Pi(a;,0 - P 2 (a:,£))] 


dx 


Applying integration by parts twice and using the boundary 
condition w^O) = 0 yields 


- - (o«(l)JV(l) + o(l)JV«(l)) «;(l) a 

+ \ (a*(0)JV(0) + a(0)JV«(0)) w( 0) 2 
+ Wa,(l)a(l)iV(l)w(l). 

Since a(x ) > a > 0 and {N, Pi, Pi} G 3d 1 ,d 2 ,e 1 ,e 2 , we have 
a(x)N(x) > Qei. Thus, by application of Lemma [5] we get 

— [ w x (x) 2 a(x)N(x)dx < ——r aei [ (w(x) — w(0)) 2 dx. 

Jo 4 J 0 

Therefore, we conclude that 

7T 2 f 1 

+ —aeitu(O) / w(x)dx 

J Jo 

-±(a x (l)N(l)+a(l)N x (l))w(l) 2 

+ ^ ^as(0).lV(0) + a(0)N x (0) — ^-aei^ w(0 ) 2 
+ Wa:(l)a(l)JV(l)«;(l). (46) 

Similarly, applying integration by parts once 


^ Jo Jo (JJ? [o(*)A(®,0]) w{£)dtdx 

\ J o w ( x ) J (J^[ a ( x )P2( x ,£)]^w(£)d£dx 

^ Jo Jo (Jx 2 w (^ d Z dx 

- J w{x) J (JP[a(x)P 2 {x,£)]\w(£)d£dx 

-w( 1) [ (d x (l)Pi(l,x) + a(l)Pi iX (l y x))w(x)dx 
Jo 

+ w(Q) j (a x (0)P2(0, x) + a( 0 )P 2 ,x( 0 , x)) w(x)dx 
Jo 

+ w x (l) / a(l)Pi(l,x)w(x)dx. 

Jo 

Changing the order of integration in the last two double integrals, 
switching the variables x and £ and using Pi(x,£) = P 2 (£,x), 


T 3 = 


l a ( x ){Pi(x,£) - P 2 (iE,f))] 


f== 


dx 


ld^_ 

2 dx 2 
1 d 2 


[a(a:)Pi(*,0] w{£)d£dx 


/ w(x) 

Jo 

n x 

w(x) 

. 

+ Jo Jo W ^[^^^ a ^ Pl ^ X ’^) W ^ d ^ dx 
+ Jo J (JJx 2 W (Q d t dx 

+ Jo J {Jw [ a (£) Pz (*’£)]) W (O d t dx 

-w(l) [ (a x (l)Pi(l,x) + a(l)Pi tX (l,x))w(x)dx 

Jo 

+ io(0) f (a x (0)P 2 (0,x) + a(0)P 2}X (0, x))w(x)dx 
Jo 

+ Wi(l) / a(l)Pi(l,a:)w(a;)dx. 

Jo 


(47) 


(48) 


Applying integration by parts once and following the same proce¬ 
dure as for T 3 , 


Now, note that for (N, Pi, P 2 ) G E t j 1 ,d 2 ,ei,e 2I we have Pi(x,£) = 
P 2 (£, x) and thus Pi(x,x) = P 2 (x,x). Exploiting this property 
and using the boundary condition, we may apply integration by 
parts twice and use w x (0) = 0 to obtain 


pi px 

/ / w(x) 1 

Jo Jo 


) w{£,)d£dx 

pi px 

/ / w(x) 1 

Jo Jo 


w^d^dx 

Jo J W( ' X ' > ^ 


) w{£,)d£dx 

Jo J W ^ ^ 


w(^)d^dx 


+ u>(l) / 6 (l)Pi(l, x)w(x)dx 

Jo 

— w(0) / b(0)P 2 (0,x)w(x)dx. 
Jo 


(49) 


Finally, employing a change of order of integration as done for T 3 
and T 4 produces 


r 5 = 


f w{x) 2 N(x)c{x)dx 
Jo 

+ Jo Jo w (Od£ dx 








J w(x) Q [c(a;) + c(0] ft (*>£)) w(£)d£dx. 


(50) 


Substituting ( I46t -( l50b into ( 145b and using Definition |2] gives us 

(Aw, Sw) + (SAw, w) 

nl rx 

< / w(x) 2 Qo(x)dx + / / w(x)Qi(x,£)w(t;)d!;dx 

Jo Jo Jo 

+ / w(x)Q2(x,£)w(£)d4dx 

Jo Jx 

+ w(0 ) ^Qsw(O) + I Q 4 (x)w(x)dx 
+ w(l) ^Q 5 w(1 ) + /: Qe(x)w(x)dx 

+ ^(1) ^Q7iu(l) + /; Qs(x)w(x)dx ) . 

Finally, using the definition of operator Q, 

(Aw, Sw) + (SAw, w) 

< (w, Qw) + w(0) 3 tu( 0 ) + J Q 4 ,(x)w(x)dx 

+ w(l) ^Qsw(l) + J Qe(x)w(x)dx 
+ w x (l) ^ Q 7 w(1 ) + r Qs(x)w(x)dx 
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